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ABSTRACT 


Control systems affected by mechanical or structural 
resonance contain at least one pair of complex poles in their 
corresponding transfer functions. The effect of the complex 
poles is to introduce resonance peaks which may cause 
instability in the closed-loop system. Such systems may be 
Stabilized by the use of cascade compensators containing one 
pair of complex zeros. 

The use of a parameter plane analysis, using a function 
of the location of the complex zeros as the two variables, in 


determining the correct location of the zeros of the compen- 


response ir. the compensated system is illustrated. 

As a result of the investigation of a number of systems 
affected by resonance phenomena using parameter plane tech- 
nigues, design criteria for cascade complex zero compensators, 
based upon the interpretation of the ‘stable region' of tne 


@oemeeometer plane, are presented. 
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ty. LNERODUCELOh 


Structural resonances in mechanical systems which, for 
example, may be caused by the twisting of motor or load shafts 
can be modelled by the introduction of at least one pair of 
complex conjugate poles in the transfer function of the system 
model. 

The effect of the complex poles is to introduce a resonant 
peak whose magnitude is dependent upon the damping ratio, Cf. 
The general form of an open loop transfer function containing 


a pair of complex poles may be written in the Bode form as: 





SES (eae et ee ere (1-1) 
(jw) Io (jwt,+1l)j{——) + 3 + 
_ K 4) Ww 
=0 el n 


The magnitude of the pair of complex poles is: 





5 A = 2 eee (1-2) 
Ge) + 3B) a] Aa wn hE 
n n gh 








whose maximum occurs where the ration of W/O, is found to be 


w/w, = V1 - 2e- doe 2 8 10g Ges) 


and the 


etek segue ry (1-4) 


arV/1-c2 


Me resonance or peak value occurs for ¢ > 0.707. 
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A plot of the peak magnitude expressed in dB is shown 
in Fig. 1.1 which illustrates that resonance peaks with 


magnitudes as high as 50 dB may be encountered. 
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Figure 1.1 


Peak magnitude as a function of damping ratio, f 


If the system gain has to be maintained at a relatively 
high level in order to satisfy other criteria, the inherent 
resonance peaks of the complex poles may cause instability 
if the resonant frequencies are contained in the bandwidth 
of the system. 

Stability may be ensured by the introduction of a cascade 
compensator whose effect is to cancel or reduce the resonant 
peaks produced by the complex poles. A very effective 


technique, if physically realizable, is to replace undesirable 
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poles with desirable ones by the introduction of a cascade 
compensator. Since the undesirable poles under present 
consideration are complex, an ideal solution would be obtained 
1f the complex poles are cancelled exactly by complex zeros 
introduced by the compensator. Unfortunately, in practice, 

it is usually physically impossible to build this type of 
compensator because the exact location of the complex poles 
can not be predicted and in fact the complex poles may change 
position during the actual operation of the system. 

The location of the complex zeros must, therefore, be 
chosen in a manner which will achieve: 

1. stability of the system, 
2. desired system response, and 
3. be physicailiy and economicaiiy realizabie. 

Lamdazuri [1] showed that a bounded region on the 'S' 
plane can be constructed with the property that, if the 
complex zeros are located within this region, stability can 
be assured. He arrived at this stability region by a method 
G@etrial and error using digital computer simulation. By 
further computer simulations of selected locations of the 
complex zeros within the stability region, he arrived at 
@ertain design criteria. 

The intent of this thesis is to arrive at this stability 
region using parameter plane techniques and to obtain design 
Gmleceria directly from the parameter plane data thus 
alleviating the need for trial and error simulation other 


than as a final check on specific system performance. 





Ti. THE PARAMETER PLANE 


The values of the parameters of a dynamic system determine 
its response to various stimuli Since the parameter values 
specify the roots of the characteristic equation of the 
a octeM. 

The first parameter plane diagram waS constructed by 
Vishnegradsky showing the dynamic operating characteristics 
as functions of the coefficients of the characteristic 
Semaelon. In 1959, Mitrovic [2] introduced a method for 
determining the roots of the characteristic equation when the 
coefficient:s of the two lowest ordered terms vary indepen- 
Siljak [3,4] to provide present day parameter plane theory. 

In general, varameter plane methods establish a relation-~ 
ship between the values of two chosen parameters, a and 8, 
and the characteristic equation of a linear dynamic system 
Beeeeransforming the characteristic Seen and using the 
fact that the transform variable, s, is complex, requiring 
that the real and imaginary parts of the characteristic 
equation go to zero independently. Two equations are thus 
determined whose simultaneous solution provides values for 
the parameters a and 8 expressed as functions of the roots 
ot the characteristic equation. 

A number of digital computer programs have been developed 


which solve the simultaneous equations and provide a graphical 


I 





relationship between the parameters and the roots of the 
characteristic equation in the form of curves on the a and 

8 plane. When chosen 's' plane contours are mapped onto the 
a, 8 plane a variety of analysis and design problems can be 
solved as a function of the two specified parameters. 

A most immediate and useful application of the parameter 
plane method is the study of stability. By mapping the 
imaginary axis of the 's' plane onto the parameter plane, 
stability regions can be defined. It is this property which 
will be used extensively in the investigation to follow. 

Although the parameter plane displays the relationship 
between the roots of the characteristic equation and the two 


parameters of the system from which an acceptable root 


pattern may be selected, it does not indicate other performance 


characteristics of the system. To ensure that at least one 
performance specification is satisfied for all values of a 


and 8, it may be imbedded in the characteristic equation by 


direct substitution and algebraic manipulation. The parameter 


plane based on the imbedded form of the characteristic 
equation will then guarantee that the desired performance 
specification is met for all values of a and 8. The desired 
performance specification adopted in this investigation is 
that the error coefficient remains unchanged by the introduc— 
tion of a cascade compensator. Hence this specification will 
be imbedded in the characteristic equation prior to obtaining 


parameter plane curves. 


vs 





Ti TI. CASCADE’ COMPEEX SZ ee) Cover SATORS 

A. 'S' PLANE "STABILITY CIRCLES" FROM PARAMETER PLANE 

DATA 

A system subject to mechanical resonances may become 
unstable if the system gain has to be maintained at a rela- 
tively high level. lLandazuri [1] has shown that cascade 
compensators may be used to stabilize the system. Ideal 
compensation is achieved when the complex zeros of the 
compensator are superimposed on the complex poles of the plant 
resulting in exact cancellation of the effect of the nesenaneee 
produced by the complex poles. Adequate compensation can, 
however, be achieved if the complex zeros are located ina 
position other than coincident with the complex poles. 

Landazuri arrived at a region of the 's' plane bounded 
byeecircles of stability” with the property that the location 
of complex zeros within this region guaranteed stability of 
the system. However, a rather difficult procedure of trial 
and error using root locus methods for specific complex zero 
locations was used to determine the limit of stability at a 
specified loop gain. The use of parameter plane methods 
simplifies the procedure considerably. 

To establish a basis for comparison, the basic system 
used by Landazuri will be reexamined. A block diagram 


representation of the system is given in Fig. 3.1, where 
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Compensator Plant 





Sas Se 
Gc (S) G(s) 
Figure 3.1 
Block Diagram of Basic System 
the transfer function of the plant and compensator is 
defined as: 
G(s) ee (3-45 
Sse t Ont So 1) 
and 
2 2 


+ ) S + 
Ss 20, w 


SEME) Sear s) ee 


(2-2) 


The real poles of the compensator were chosen such that 
their contribution to the transient response of the overall 
system will disappear before the maximum overshoot is reached. 

The value of the gain, k, at the limit of stability was 
calculated to be 0.1. For this value of gain the error 


@eetficient, Kv, is 0.l. 


i> 





The system transfer function is 


x (s* a 25, Ss + w 4) 
GcG(s) = ———_-————- SSS SS (3273) 


s(s+10) (stl? 2 Onlmeme 1) 





and the error coefficient, Kv, is 





x n 
Ken = 500 : (3-4) 


If the error coefficient is imbedded, equation (3-3) becomes 


200Kv/u,* 5 PERO c/o, S + 200Kv 
nae ane a ee 
s(st+t10) (s+20) (s°+0.1s + 1) 


and defining a = L/w,” and B= C/w,, 
allows equation (3-5) to be written as 


200Ky a s* + 2(200Kv) 8 s + 200Kv 
GeG(s) = ———_________, —__-__ (oe) 
s(s+10) (s+20)(s~ + 0.1s + 1) 





Mie characteristic equation of the system: 
1 + GceG(s) = 0 (2=7)) 


will, therefore, contain the two parameters a and 8. 

The results of implementing equation (3-7) ae a parameter 
plane problem for values of Kv 10, 50 and 100 percent above 
its value at the limit of stability of the uncompensated 
System are shownin Fig. 3.2, 3.3 and 3.4 respectively. 

A mapping of the imaginary axis of the ‘'s' plane is 


represented by the zeta = 0.0 curve on the parameter plane, 
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1.@€., values of a and 6 chosen to lie on the _ = 0.0 

curve ensure that cne pair of complex roots of the system 

is located on the imaginary axis of the 's' plane. The system 
is, therefore, at the limit of stability if all of the remain- 
Mey roots are in the left-half plane. From Figs. 3.2, 3.3 

and 3.4 values of a and @ that lie on the respective zeta = 0.0 
curves define, therefore, the stability boundary on the 's' 


plane in terms of the OO, and ¢ of the complex zeros since 


tees dey shots, (6 = | itt : (3-8) 


The resultant stability boundaries are shown in Fig. 3.5 
for three values of Kv. It should be noted that these 
boundaries are identical to the ones determined bv Landazuri 
with the exception that they should not be referred to as 
Stability circles since their actual shape deviates consider- 
ably from circular as illustrated particularly by the results 
obtained from Kv = 0.20. 

The use of parameter plane date has considerably simpli- 
fied the procedure of obtaining 's' plane stability bound- 
aries. A considerable amount of labor remains, however, in 
converting parameter plane data for use on the 's' plane. 
Mies 15, in fact, unnecessary since the parameter plane 
curves can be drawn in terms of the on and © of the conpaex 
zeros directly. The resultant graphs impart as much, and in 


fact more, information as the 's' plane without intermediate 


steps as will be shown in the following section. 
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meme: a0 = 1.0/fin, 6 = 0.5/in 


Figure 3.2 


Parameter Plane Curves for a = 1/w 


a and € = c/w, 








meee: © — 1.0/in, 6 = 0.5/in 
Figure 3.3 


Parameter Plane Curves for a = L/w, and 6 = C/u, 
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is 
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Beale: a9 = 1.0/fin, 8 = 0.5/in 


Figure 3.4 
Parameter Plane Curves for a = L/w? and B = C/O, 
Kv = 0.20 
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B. STABILITY REGIONS ON THE PARAMETER PLANE 

The stability boundaries as shown in Fig. 3.5 clearly 
define that region on the 's' plane in which a pair of complex 
zeros may be located to guarantee stability of the system. 
Since these boundaries were determined from parameter plane 
data, the parameter plane itself can be manipulated to show 
the stability region just as clearly as the 's' plane. 

The two parameters a and 8 were previously defined as 
functions of the location of the complex zeros, namely L/w," 
and C/o, respectively. This choice was imposed by the fact 
that the error coefficient was imbedded in the transfer 
function of the system. Parameter olane curve calculations 
must, ‘therefore, be based upon these two parameters. For 
MeGiecang Purposes, however, any function of the original a 


and 8 may be used. One obvious choice was to define two new 


macbemeters 
A f- A 
a, = V1/oa and Sq BY1/a (3-9) 
in terms of a and 8 with the result that ai Sites and By = C 


of the complex zeros. 

Figures 3.6, 3.7 and 3.8 were plotted with these new axis, 
mueeotherwise are identical to Figs. 3.2, 3.3 and 3.4 respec- 
tively. The stability region is clearly defined by the 
zeta = 0 curve indicating that values of a and 8 chosen to 
lie on this curve will guarantee that one pair of roots of 
the closed-loop system will lie on the imaginary axis of the 


's' plane and, if the zeta = 0 curve forms an enclosed 


Ze 





region on the parameter plane, that the remaining roots are 
meeene left half of *ther 4s “Wallane. 

Constant zeta curves of values other than zero can, of 
course, also be plotted as shown. This means that values of 
a and 8 chosen to lie on these constant zeta curves will 
guarantee that one pair of roots of the closed-loop system 
will have a corresponding zeta and, if the constant zeta 
curve lies within the stability region, that the system is 
stable; i.e., all roots are in the left half of the ‘'s' plane. 
The constant zeta curves can, therefore, be used as a tool 
in meeting design specifications since they define values of 
a and 8 which will Posuake in a specified zeta. Unfortunately 
they do not guarantee dominance of the selected roots and 
Shovld, therefore, only be used as a rough guide and net as 
an absolute indication of the closed-loop system performances. 
The problem of dominance of selected roots will be further 
discussed in sections to follow. 

To illustrate, however, the information available from 
the constant zeta curves, two sets of veoce OL Wn and t were 
chosen as indicated by points A and B on Fig. 3.7. Clearly 
the system response for values of Ww and ~ defined by point 
B should be more damped than the response for values of w 
emer c defined by point A. That this is, in fact, the case 
fomiliustrated in Fig. 3.9 and 3.10 which represent the 


transient response of the system to a step input for values 


fore Ow and z defined by points A and B respectively. 





Since the stability region exists for values of |¢| > 1 
as well, it. is immediately obvious that the system can be 
eat lized mot only by the addition Of a pair of complex 
zeros, but also by two real zeros, which was not obvious from 


the 's’ plane stability region. 


For Kv = 0.15, values of vo, = 0.5 and ¢ = 4.0 were chosen 
which result in the location of two real zeros at s = -0.063 
and s = -3.94 as shown in Fig. 3.12. That the resultant 


system is stable is immediately obvious from Fig. 3.12 and 
illustrated in Fig. 3.11 which represents the transient 
meoponse of the system for the chosen values of WO and Cf. 

Since it 1S considerably simplex and more economical to 
build real zero compensators than complex ones, the practical 
designer would certainly prefer their use if possible. The 
parameter plane, therefore, allows nim to determine whether 
real zero compensators can be used for a particular system. 

Landazuri [1] arrived at design criteria which stated 
that the complex zeros should be located so that: 


mc. = > , and 


The parameter plane curves do not directly indicate criteria 
Mee but criteria 2. is obvious from the curves. Since the 
values of wo for curves of Constant 2ora 7 02) areud eiese 
than 1.0 which is the value of oe a value of w < 9 will 
result in a more damped transient response. 

It has been shown that the parameter plane can display 


stability regions defined by the zeta = 0 curve. The 
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Figure 3.6 


Parameter Plane Curves for a 
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Figure 3.9 


Transient Response for Point A. 
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Transient Response for Point B. 
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Figure 3.11 


Transient Response for Compensated System with Real Zeros. 
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parameter plane canwalso give jqurdancem, 16cm the 
Semplex zeros to ensure a Specvitvedseype of respemico by 
using the constant zeta curves. It would appear, therefore, 
that the parameter plane is ideally suited as an aide to the 
designer faced with the problem of complex zero compensation. 
So far, however, only a very simple system has been investi- 
gated. The following section will, therefore, generalize 
the problem of cascade complex zero compensators and investi- 
gate a number of more complicated cases such as: 

1. systems of higher order, 

2. effect of round-off poles, 

3. systems with two resonance peaks, and 


4, systems with resonance-antiresonance doublets. 


C. PARAMETER FPLANE DESIGN OF COMPLEX ZERO COMPENSATORS 
1. The General Equation 
The block diagram representation of the general 
system under consideration is illustrated in Fig. 3.1. The 
generalized transfer function of the system may be expressed 


in the following form 


Je m 5 5 
; 1 please War 
k (s*+2tu roan 2) K ee +2n sus us ] 
GcG(s) = —_—_——"— Fei ea eee ene SS 
= it 2 2 
ae! sN I (stp, ) I [s +2n, up ty, J 
k=0 k=0 
(SILO) 
where: a. K is large enough to cause instability in the 


uncompensated system, 


oe Zor Pye MT eis ib and Wy, are known, 


BZ 





c. p is chosen sufficiently large so that its 
contribution to the transient response of the 
compensated system is negligible, 


2 
qd. k= tas so that the error coefficient, Kv, 


remains unchanged, and 
e. mn and uw are equivalent to the generally accepted 
interpretation of t and a, respectively and are 
used only to differentiate the zeros of the 
compensator. 
In order to transform equation (3-10) into an accep- 
table parameter plane form the two parameters a and 8B are 


Gefined as: 


A 2 
Q = 1/0, and B 


and upon substitution the general transfer function becomes: 


», We ss mn © 2 
Ione s +267stl) i (stz,) qT (s +2n.UstHs”) 
Gee(s) = ake ae (3-11) 
Nee a ree 2 
s I (stp, ) (stp) I (s F2n Hy tH, ) 
k=0 k=0 


meam which the characteristic equation (3-7) can be formed: 

To illustrate theuse of the parameter plane in the design 
of a suitable complex zero. compensator for a system affected 
by mechanical resonances, a number of examples will be studied 
which contain various combinations of pole and zero terms 


included in equation (3-1l). 
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The mathematical solution for the parameter plane 
Curves will, of course, allow values of @ and @ which may be 
both positive or negative. Physically, both positive and 
negative values of 8 are permissible since a positive 8 
indicates a complex pair or two real zeros in the left half 
of the 's' plane, while a negative 8 indicates zeros in the 
right-half plane. Negative values of a, however, are not 
permitted since no corresponding physical interpretation 
exists. Negative values of a are, therefore, ignored in the 
interpretation of stability regions. 

Whenever negative values of oa occur in the parameter 
plane solution, the previously defined transformation (3-9) 
presents obvious difficulties. The problem can, of course, 
Bemocsolved by setting a, = 0 whenever a is negative or by 
taking the square root of the absolute value of l/a. Digital 
computer plotting of parameter plane curves under these 
conditions introduces extraneous lines whenever a becomes 
negative so that the clarity of the parameter plane presenta- 
tion diminishes considerably. For this reason the following 
parameter plane curves are plotted with a = L/w,” and g = C/u), 
without further transformation. 

The stability regions are still as clearly defined 
as with the previously defined transformation, but negative 
welues of a must be ignored, and it is not immediately 
obvious whether a chosen combination of a and 8 will result 
in a pair of complex zeros or in two real zeros. The advantage 


of clearly defined regions of stability obtained by plotting 
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a and 8 directly outweighs the indicated iadvaintadee and, 
hence, in tne following parameter plane analysis a and 8 will 
be plotted without further transformation. 
Za onalysis of a FOuUrEnSerdermoy seems 
The type 1 - 4th order system under investigation 
has an uncompensated transfer function 


CB) a (3-2, 


s(st2) (s“+0.1 s+1) 
aueeats root locus is shown in Fig. 3.13. 

The value of the gain, K, at the limit of stability 
was found to be 0.236 and the corresponding error coefficient, 
Kv, is 0.118. The error coefficient was increased by 50 
percent to a vaiuc cf 9.177 and the system was stabilized by 
the introduction of a cascade complex zero compensator whose 


easter function is: 


100/w,* (s*+270, a wi) 
Ge(s) = ——~—__, —*_"— (3-13) 
(s+10) 
The resultant parameter plane curves are shown in 
moe 3.14 with the stability region clearly defined as the 
area bounded by the zeta = 0 curve. Bounded regions are also 
@emined for zetas equal to 0.1, 0.2 and 0.5 which are entirely 
contained within the stability area. Thus values of a and 8 
chosen to lie within these bounded regions would guarantee 
that all the roots of the closed-loop system would have a 


corresponding zeta or higher. 
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From Ene parameter pilane aavems alles of o = 2.5 and 
@ = -1.15, which lie within the area bounded by zeta = 0.5, 
should give the system a transient response which is most 
acceptable from the standpoint of maximum overshoot and 
settling time. That this is, in fact, the case is shown in 
Fig. 3.15 which represents the transient response of the 
System for the indicated values of a and @. The corresponding 
meee Locus is shown in Fig. 3.16. 

Two other values of a and 8 were chosen to establish 
a basis of comparison. The transient response for a = 1.0 
and 8 = 9.0 is shown in Fig. 3.17 and the corresponding root 
Mecus in Fig. 3.18. From Fig. 3.18, it is clear that almost 
exact cancellation of the complex poles has taken place which 
accounts for the highly damped response indicated in Fig. 3.17 
for the chosen value of gain. This type of response was not 
directly obvious from the parameter plane and is due to the 
fact that almost exact cancellation has taken place with the 
resultant residues due to the complex poles being extremely 
small. 

The oscillatory response as shown in Fig. 3.19 for 
a = 16.0 and g = 0.0 was expected from the parameter plane 
data although its exact nature was not apparent since it is 
dependent upon the residue of the complex poles at the 
Bereticular gain. 

It was previously mentioned that, unfortunately, the 
parameter plane does not give an indication of dominance of 


the closed-loop system roots selected by a particular choice 
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Figure 3.13 


Root Locus of Uncompensated 4th-order System. 
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Figure 3.15 


Transient Response of Compensated 4th-order System 
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Figure 3.16 
Root Locus of Compensated 4th-order System 
Qo = Ze"; B = -1.15 
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Figure 3.18 
Root Locus of Compensated 4th-order System 
o_o 6 = O20 
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Figure 3.20 
Root Locus of Compensated 4th-order 
a = 16.0, 8 = 0.0 
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of a and 8B. Such an indication an only come from the root 
Hocus of the compensated Syetem takingeintoudcecoun men] 
selected gain. If the transient response obtained for a 
particular choice of a and 8 is not suitable, the root locus 
of the corresponding system will give an indication of why 
this should be the case. Relocation of the zeros on the root 
locus may lead to a more Suitable response and the selected 
zero location may be checked on the parameter plane to 
determine that stability of the system is still ensured. It 
may, therefore, be necessary to use both the parameter plane 
and the root locus of the system tc predict a suitable 
response. The parameter plane alone, however, will still 
indicate the location of compensating zeros to give, at least, 
a first approximation to a desired response. 
wae Nc Wse OF Keung orf Poles 

Mrepractieal engineering the use Of Ene) round or: 
pole is widely accepted as a means of stabilizing a system 
which is basically unstable due to mechanical resonance 
peaks. The effect of a round-off pole can best be illustrated 
on a Bode diagram such as shown in Fig. 3.21. If an addition- 
al pole is introduced into the system at a frequency less than 
the undamped natural frequency of the complex poles, the 
value of the magnitude of the resonant peak is decreased. 
If the reduction in magnitude is sufficient to ensure that the 
resonant peak does not cross the zero dB axis, the stability 


of the system will be assured. 
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Bugure 3.21 


Bode Diagram Representation of the effect of Round-off Poles 


To illustrate the effect of a round-off pole, an 
additional pole at s = -0.5 was introduced into the 4th-order 
system described previously. The uncompensated system 


transfer function becomes, therefore, 


G(s) = ——_____+_,____ Ger) 
s(st0.5) (s+2) (s°~+0.1 stl) 


mies value of K at the limit of stability has been increased 
to 0.619 with a corresponding increase in the value of the 


pmmor coefficient, Kv, at the limit of stability to 0.619. 
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When compared with the value of the error coefficient ene aleve. 
fea-order system at the limptyeces cabs lay sane nem eesc met 

525 percent is allowed before the system describec by equation 
3-14 becomes unstable. This is illustrated in Fig. 3.22, the 
meet locus of the system. Sinee the reot loci emerging from 
Ehe complex poles curve towards the origin before crossing 

the imaginary axis a larger value of gain may be applied to 
the system before the limit of stability is reached. For a 
value of Kv equal to Q.177, as in the previous example, the 
system is stable as shown in Fig. 3.23. 

The introduction of a complex zero compensator of the 
form described by equation 3-13 is not required to stabilize 
the system Since it is already stable, but may be used to 
modify the system so that a more suitabieé transient response 
can be achieved. The parameter plane curves for this system 
fea constant value of oreeaa1 to 0.177 are shewniein Ber 
3.24. The stable region is defined by the zeta = 0 curve for 
positive values of a. The parameter plane curves immediate-y 
show that complex zeros would have to be introduced Dri Maire, 
in the right half of the 's' plane in order to bring the 
Syotem to the limit of stabilitv at the given value of gain. 

An inspection of the curves of Fig. 3.24 would lead 
to the conclusion that the optimum response of the system 
could be expected for values of a and 8 in the neighborhood 
of 1.0 and -1.5 respectively, which incidentally give rise to 
two real zeros. The transient response for these values of 
o and 8 is shown in Fig. 3.25 and the corresponding root locus 


in Fig. 3.26. A comparison of Figs. 3.23 and 3.25 indicates 
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Root Locus of Uncompensated System with Round-off Pole 
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Figure 3.24 


Parameter Plane Curves of System with Round-off Pole 
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Frgure» 3.29 


Transient Response of Compensated System with Round-off 
Pole a= 1.0, 8 = -1.5 
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Figure 3.26 
Root Locus of Compensated System with Round-off Pole 
fe= 1.0, 8B = -1.5 
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that no Significant advantage is gained in the transient 
response by introduction of the complex zero compensator, in 
fact, both maximum overshoot and settling time have been 
increased. Exact, or at least nearly exact, cancellation of 
the complex roots by compiex zeros does, however, improve 
mies transient response as shown in Figs. 3.27 and 3.28, 

The use of round-off poles alone, may be preferrable 
to the use of complex zeros in stabilizing a system affected 
by resonance peaks because it is a much more economical 
method of achieving the desired results. Unfortunately the 
method is jiimited to systems that have a sufficiently large 
phase margin at the saben cross-over, in terms of the Bode 
meaeram, SO that the reduction in phase margin introduced 
by the round-off nole does not drive the basic s 
No particular advantage is achieved by the use of a combin- 
ation of round-off poles and complex zero compensators. 

Exact, or nearly exact, cancellation of tne complex poles may, 
of course, be used, but this would have been just as effective 
without the introduction of the POae ore eels 

4. Analysis of a System with two Resonance Peaks 

In practical control systems containing mechanical 
linkages or shafts, two or more resonant frequencies are often 
encountered which may be medelled by the introduction of the 
appropriate number of complex poles in the transfer function 
of the system. The feasibility of compensating a system with 
two resonant peaks by the use of one pair of complex zeros 1s 


investigated in the following paragraphs. 
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Théestranstfer Lunctiommet Ene wsyorem containing two 


pairs of complex poles is given by: 


G(s) = 


K (s+2) 


Bist 5) (oto) (oneal Soul) ee 


(3-15) 
+0.3 st4) 


the root locus of the system described by equation 3-15 is 


Sewn in Fig. 3. 


Whee Er 


the velwesor Keac Che. limit Vor Ss talneiney 


was found to be 1.11 and corresponds to a value of Kv equal 


Tom. 222. 


Adding a complex zero compensator of the form 


described by equation 3-13, and increasing the value of Kv 


by 50 percent to 0.333 results in the parameter plane analysis 


Siew in Fig. 3.30. 


is bounded 

LO 
and 8 must 
Bena first 


the region 


transient response which 1s most acceptable. 


by the zeta = 


be chosen which lie within 
approximation a value of a 


encllosed by the zeta = 0.2 


As previously explained the stable region 


O curve for positive values of a. 


Oi S Gen, ae ee oe 
the stability region. 
and 6 which falls within 


cumver should provaide a 


The transient 


response for values of a and 6 equal to 0.07 and -0.7 respec- 


tively is shown in Fig. 


ig. 3. 


325 


It is obvious from Fig. 


3.31 and the associated root locus 


Seedy that altnougm 


the settling time may be acceptable, the maximum overshoot 


would appear to be excessive. 


For values of a and 8 equal to 


0.2 and -0.2 respectively the transient response is shown in 


Pare 3-33 and the corresponding root locus in Fig. 3.34. 


AS 


expected for these values of a and 8, the system is less 


damped so that the settling time has been increased, but the 


maximum overshoot has been decreased to an acceptable level. 
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Clearly a compromise between the two chosen values of o and 
8B is indicated to obtain the optimum response under the given 
conditions. 

Since stability was achieved by the use of two real 
zeros rather than a pair of complex ones for the choice of 
a and 8 as shown in Fig. 3.32, it may be economically advan- 
tageous in this case to chose values of a and 6 which will 
result in two real zeros. The parameter plane will, of 
course, indicate what range of values of a and 8 is permissible. 

It should also be pointed out that cancellation of 
one pair of complex poles is no longer feasible and, in fact, 
prohibited since an unstable system will result from such an 
action. The parameter plane curves indicate this clearly 
Since only values of a between 0.95 and 0.22 for 8 = 0, which 
correspond to Wn between 4.47 and 2.13 for t = 0, are within 
the stable region. 

A system with only one resonance peak can, of course, 
always be stabilized by exact or nearly exact cancellation of 
the complex pole and the parameter plane analysis may not be 
required. For systems with two or more resonance peaks, 
however, cancellation of one pair of complex poles with one 
pair of complex zeros may not ensure stability, but a para- 
meter plane analysis will immediately indicate if stability 
can be assured with only one pair of complex zeros and, if 


so, Where they should be located. 
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Root Locus of Uncompensated System with Two Resonance Peaks 


58 





gaa 





Seale: of = 0.5/fin, 8 = 0.3/in 


Figure 3.30 


Parameter Plane Curves for System with Two Resonance Peaks 
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Transient Response of Compensated System with Two Resonance 
Peaks a= 0.07, 8 = -0.7 
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Root Locus of Compensated System with Two Resonance Peaks 
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Figure 3.34 
Root Locus of Compensated System with Two Resonance Peaks 
oes 0.2, 68 = -0.2 
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5. Analysis of Systems with Resonance-Antiresonance 
Doublets 


Although the effect of structural resonances in 
practical control systems may be modelled by the introduction 
of a number of complex poles in the transfer function of the 
system as described in the previous section, a more accurate 
description of the system is often found by modelling the 
structural resonances as resonant-antiresonant doublets. The 
term 'doublet' 1s used to indicate that the characteristic of 
the structural resonance is such that a resonance peak 
combined with an anti-resonance trough appear at closely 
related frequencies in the frequency response of the overali 
system. Such a doublet can be modelled by the inclusion of a 
combination of one pair of complex poles and one pair of 
complex zeros in tne system transfer function. 

To illustrate the effect of structural resonances 
which can be modelled in the form cf resonance-antiresonance 
Gemet@ets, a basic type 1 - 5th order system whose transfer 


function is, 


K(s+0.5) 


G(s) = 5 
s(s+0.2) (st+4) (s~+4s4+13) 


(3= Gy 


and whose corresponding root locus is shown in Pig@assss> 

will be subjected to a number of structural resonances in the 
form of doublets. For a selected gain, K, equal to 31.2, the 
error coefficient, Kv, is 1.5 and the system transient 


response, as shown in Fig. 3.36, was deemed acceptable from a 
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standpoint of maximum overshoot and settling ne An 
analysis of this system when subjected to structural resonances 
will be presented in the following sections. 

a. A System with One Doublet 

The inclusion of one resonance~antiresonance 

doublet in the system described in the previous paragraph will 
result in the addition of one pair of complex zeros and one 
pair of complex poles to the transfer function of the system 
which, therefore, becomes: 


Z 5 
mie) = K(st0.5) (s~+0.12s+1.22) (3-17) 


s(s+0.2) (st4) (s7+4s+13) (s“+0.1 s+0.8) 
@iemroot locus of this system is shown in Fig. 3.37, indicating 
f@etwfor the specified value of Kv the system us unstable. In 
order to stabilize this system a complex zero compensator of 
the form described by equation 3-13 was used and a parameter 
plane analysis for the specified value of Kv was carried out. 
The resultant parameter plane curves are shown in Fig. 3.38 
indicating that a stable region exists, as shown by the shaded 
area. 

Since a portion of the area bounded by the zeta = 
0.2 curve lies within the stable region, a selection of a and 
8B from that area of intersection should result in a compensated 
system which has the most acceptable transient response achiev- 
able under the specified conditions. For values of a and B 
equal to 1.0 and 0.0 respectively, which are contained within 


the area of intersection, the transient response of the 
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Parameter Plane Curves for System with one Doublet 
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Transient Response of Compensated System with one Doublet 
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70 








Sete: xX and y = 2.0/in 


Figure 3.40 
Root Locus of Compensated Systems with one Doublet 
a= 1.0, 8 = 0.0 
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Figure 3.41 
Transient Response of Compensated System with one Doublet 
a= 2.0, 8B = - 0.15 
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Figure 3.42 
Root Locus of Compensated System with one Doublet 
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Figure 3.43 
Transient Response of Compensated System with one Doublet 
Oe= 92.0, 6 = 1.0 
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compensated system is shown in Fig. 3.39 and the corresponding 
root locus in Fig. 3.40. A comparison of the transient 
mesponse shown in Fag. 3.39 with those displayed ian verge. 
3.41 and 3.43 would support the assumption that the most 
acceptable transient response is obtained by choosing values 
of a and 8 which lie within the area of intersection of the 
stable region and the area bounded hy the largest constant 
zeta curve. The transient response shown in Fig. 3.39 appears 
acceptable as far as maximum overshoot is concerned, but 
displays a certain amount of small amplitude ‘'ringing' which 
may or may not be acceptable in a particular system. The 
ringing is not present in the transient response shown in 
Fig. 3.41 and, thus, a compromise between the two values of 
a and 8 selected may be indicated to obtain the optimum 
response of this particular system. The parameter plane 
curves, however, do indicate what values of a and 8 will 
Stabilize the system and, to a first approximation, what 
values will give the most acceptable transient response. 

b. A System with Two Doublets 

The inclusion of a second resonance - antiresonance 

doublet in the system described previously gives rise to a 


system transfer function of the form: 


K(st0.5) (s°+0.12s¢1.22) (s“+0. 35+9.0) 


ae10.2) (e+4) (o> 245413) (s°20.1 st0.8) (S7*+0.2S+7.85) 
(3.18) 


G(s) = 


The root locus of the system described by equation 3-18 is 
shown in Fig. 3.45 and again the system is unstable fora 


value of Kv equal to 1.5. As before, it was attempted to 
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Stabilize the system using a complex zero compensator of the 
form described by equation 3-13. The resultant parameter 
mlane analysis, shown in Fig. 3.46, again indicatesmamctapre 
region closely resembling the one obtained for one doublet. 
Since the instability in this particular system is caused by 
Ene first doublet, as shown in Fig. 3.45, it is reasonable 
to assume that the system can be stabilized in a manner 
Similar to the one used for one doublet. The similarity of 
the two parameter plane curves, Figs. 3.38 and 3.46, substan- 
tiates this assumption. The small differences that can be 
detected between the two parameter plane curves are due to 
the presence of the extra doublet. 

As in the previous example, the most acceptable 
transient response of the compensated system should be 
obtained by selecting values of a and 8 which lie within the 
area of intersection of the stable region and the area 
enclosed by the zeta = 0.2 curve. The transient response of 
the compensated system for values of a and 8 within the area 
of intersection under consideratior. is shown in Fig. 3.47. A 
Semoadrison of Fig. 3.47 with Figs. 3.49 and 3.51, which 
represent the transient response of the system for other 
combinations of a and 8, indicates that the assumption as 
stated was valid. Again verifying that a parameter plane 
analysis not only indicates whether a system can be stabilized 
by a selected method, but also what range of values of the 
parameters will result in a most favourable response of the 


compensated system. 


17 








+ 
B52 


+ 
+ + 
Sgt 


u ne a a ee 


Se@ere: xX and y = 1.0/in 


Pargigeess os5 


Root Locus of Uncompensated System with Two Doublets 


78 





ba 


-3c6 





Seete: a = 3.0/fin, 8 = 2.0/in 


Figure 3.46 


Parameter Plane Curves for System with Two Doublets 
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Figure 3.48 
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Transient Response of Compensated System with Two Doublets 
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Figure 3.52 
Root Locus of Compensated System with Two Doublets 
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c. A System with Three Doublets 

In the two previous examples of the effect of 
resonance - antiresonance doublets the undamped natural 
frequency of the complex poles was less than that of the 
accompanying complex zeros. This is, of course, not always 
the case and for that reason the undamped natural frequencies 
of the third doublet included in the system under investiga- 
tion have been reversed, i.e., the undamped natural frequency 
of the complex zeros is less than that of the complex poles. 

The transfer function of the uncompensated system 
with three doublets becomes: 


2 
ccs) = ____K(st0.5) (s +0 .12stl.22) 0, 


Ss (st0.2) (S44) (s- te Ase sme eto ee tcIC) 


(s*+0.3s+9.0) (s°+0.88+19.5) 
(e-t0 ety 85) (ac tenosme 


and the root locus of this system is shown in Fig. 3.53. It 
should be noted that Fig. 3.53 shows dheie aine complex zeros 
Seeeme third doublet are located closer to the origin than 
the complex poles indicating that the undamped natural 
frequency of the complex zeros is less than that of the 
complex poles. 

Since the system, as before, is unstable because 
of the presence of the first doublet, as shown in Eig. 320m 
it is again reasonable to expect that a parameter plane 


analysis of the feasibility of stabilizing the system with 
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Parameter Plane Curves for System with Three Doublets 
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Figure 3.56 
Root Locus of Compensated System with Three Doublets 
a= 1.2, 8 = 0.0 
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Figure 3.59 
Transient Response of Compensated System with Three Doublets 
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Root Locus of Compensated System with Three Doublets 
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one complex zero compensator of the form — by equa- 
tion 3-13 should result in parameter plane curves resembling 
Eeese Shown in Figs. 3.38 and 3.46: his is; anetact cnc 
case aS illustrated by Fig. 3.54 which represents the para- 
meter plane curves for the compensated system with three 
@emplets. 

Similar to the procedure followed in the previous 
two cases, three sets of values of a and 8 from within the 
stable region were chosen and the corresponding transient 
responses of the compensated systems are shown in Figs. 3.55, 
3.57 and 3.59 together with their respective root loci in 
fee. 3-90, 3.58 and 3.60. AS CxypectCammele Mesemaeeopeaone 


transient response was obtained for values of a and 8 within 


area bounded by the zeta = 0.2 curve as shown in Fig. 3.55. 

d. Discussion of Results 

The analysis carried out in the last three 

sections have shown that a basic control system which is 
stable at a selected value of the error Pesemiotene can become 
unstable at the same value of Kv if the effect of structural 
resonances is taken into account. It has also been shown that 
Such an unstable system can be stabilized by the use of a 
complex zero compensator if the location of the pair of 
complex zeros or real zeros is such that the root loci of the 
compensated system are adjusted so that none of the roots of 
the closed loop system at the specified gain are located in 
the right half of the 's' plane. The problem of finding the 


correct location of the zeros of the compensator to achieve 


eh, 





stability of the system is solved by the use of a parameter 
piane analysis where the two variable parameters are selected 
to be functions of the location of the zeros of the compensa- 
tor. A parameter plane analysis not only indicates the proper 
location of the zeros of the compensator to achieve stability 
by the use of the zeta = 0.0 curve, but also gives an indica- 
tion of where the zeros should be located to assure a desired 
type of transient response by the use of constant zeta curves 
of values other than zero. 

The optimum transient response of a system 
affected by structural resonances and stabilized by the use 
of a cascade complex zero compensator may not approach the 


desired response achieved by the system without structural 


Figs. 3.39, 3.47 and 3.55 are not as good as the one shown 
in Fig. 3.36, but under the given conditions they are 
approximately the best that can be obtained from the stand- 
point of maximum overshoot and settling time. If an area of 
intersection between the stable region ana an enclosed area 
bounded by a constant zeta curve of value greater than 0.2 
had been indicated, a more favourable transient response 
could have been possible. The parameter plane curves shown 
in Figs. 3.38, 3.46 and 3.54 indicate, however, that this is 
not the case in the presented examples and that under the 
given conditions and for the selected compensator the optimum 
has been reached. The value of a parameter plane analysis 


lies, therefore, in the immediate availability of the location 
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of the zeros of the compensator to achieve stability of the 
system and in the information available from which a first 
approximation of the obtainable transient response can be 


formed. 
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IV. CONCLUSIONS 


Control systems affected by mechanical or structural 
resonances contain at least one pair of complex poles in the 
corresponding transfer functions. The presence of the complex 
poles causes BS QRERE peaks to occur in the open-loop frequency 
response of the system which may produce instability in the 
closed-loop system. The use of complex zero compensators to 
Stabilize such systems has been widely accepted particularly 
when complete cancellation of the complex pole by the complex 
zero of the compensator can be achieved. 

Exact cancellation by one pair of complex zeros is, of 
course, impossible to achieve if the frequency of the 
mechanical or structural resonance varies during the oper- 
ation of the system, or if resonances at two or more different 
frequencies occur. Such a system can, however, be stabilized 
Peeeeecascade compensator containing one pair of complex or two 
real zeros if the choice of zeros is such that they modify 
the root locus of the system in a manner which will ensure 
that all roots of the closed-loop system remain in the left 
half of the 's' plane. A parameter plane analysis, using 
functions of the location of the complex zeros as the two 
variables, provides an effective method of determining the 
correct location of the zeros of the compensator to ensure 
stability. At least one other important system criteria, 


such as constant error coefficient, may be imbedded in the 
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analysis to ensure that the solution obtained will satisfy 
mee eSpeCl Ered Criteria asm 

A parameter plane analysis may be used not only as a 
method to determine the correct location of compensating zeros 
to ensure stability, but also as an aid in the design of the 
compensator to achieve a desired system performance. The 
exact response of the compensated system can not be predicted 
from a parameter plane analysis, but an approximation based 
On second order system theory can be obtained if constant 
zeta curves of values other than zero are plotted on the 
Parameter plane. 

It is felt that a parameter analysis using cascade compiex 
zero compensators, as has been illustrated in this investi- 
gation, can provide considerabie aid to the practical designer 
faced with the problem’of a system affected by mechanical or 
aemuecural resonances. 

The examples presented in this investigation and the con- 
clusions drawn from them by no means exhaust the study of the 
applicability of parameter plane analysis in the design of 
complex zero compensators. The following questions are 
presented as suggestions for possible future investigations. 

1. What is effect of the location of the poles of 
the compensator, i.e., does their location 
seriously affect the size of the stable region 
on the parameter plane? 

2. Can a system containing several doublets be 


Stabilized with one pair Of complex Zerosa1 


og 





the root locus geometry srs (such Ehat seine oe 
loci of two doublets enter the right-half plane 
and assuming that, at sufficiently high gain, the 
roots of the lower frequency doublet reenter the 
left-half plane, thus causing instability of the 


system due to the second doublet? 
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